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Abstract 

In [7], Berthelin, Degond, Delitala and Rascle introduced a traffic flow model 
describing the formation and the dynamics of traffic jams. This model consists of a 
Pressureless Gas Dynamics system under a maximal constraint on the density and 
is derived through a singular limit of the Aw-Rascle model. In the present paper 
we propose an improvement of this model by allowing the road to be multi-lane 
piecewise. The idea is to use the maximal constraint to model the number of lanes. 
We also add in the model a parameter a which model the various speed limitations 
according to the number of lanes. We present the dynamical behaviour of clusters 
(traffic jams) and by approximation with such solutions, we obtain an existence 
result of weak solutions for any initial data. 

Key words: Traffic flow models. Constrained Pressureless Gas Dynamics, Multi- 
lane, Weak solutions. Traffic jams 
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1 Introduction 

Classical models of traffic are splitted into three main categories: particle models (or 
"car- following" models) ^15^ i3j, kinetic models [231 (211 EQ] [IE], and fluid dynamical 
models [HI HH [22l |2l ]26l [HI [16]. Obviously, these models are related; for exam- 
ple in [T], a fluid model is derived from a particle model. See also [T7]. Here, we 
are interested in the third approach, which describes the evolution of macroscopic 
variables (like density, velocity, flow) in space and time. Let us recall briefly the 
history of such models. 

The simplest fluid models of traffic are based on the single conservation law 

dtn + dxf{n) = 0, 

where n = n{t, x) is the density of vehicles and f{n) the associated flow. This model 
only assumes the conservation of the number of cars. Such models are called "first 
order" models, and the first one is due to Lighthill and Whitham [TH] and Richards 

m- 

If we take the fiux /(n) = nu with u = u{t, x) the velocity of the cars, we add 
a second equation of equilibrium related to the conservation of momentum. This 
approach starts with the Payne- Whitham model [211 I22j . 

But the analogy fiuid- vehicles is not really convincing: in fact, in the paper [12] , 



2 



Daganzo shown the Umits of this analogy, exhibiting absurdities which are imphed 
by classical second-order models, for example, vehicles going backwards. To reha- 
bilitate these models, Aw and Rascle proposed in [2] a new one which corrects the 
deficiencies pointed out by Daganzo. In particular, the density and velocity remain 
nonnegative. 

The Aw- Rascle model is given by 

f dtn + d^{nu) = 0, 

\ {dt + ud^){u + p{n)) = Q, 

or in the conservative form 

f dtn + dx{nu) = 0, 

\ (?t(n(u + p(n))) + dx{nu{u + p{n))) = 0, 

where p{n) ~ is the velocity offset, which bears analogies with the pressure in 
fluid dynamics. 

In fact, this model can be derived from a microscopic "car- following" model, as it 
has been shown in [IJ. But even the Aw- Rascle model exhibits some unphysical 
feature, namely the non-propagation of the upper bound of the density n, making 
a constraint such that n <n* impossible (where n* stands for a maximal density of 
vehicles) . 

Some constraints models have been developed these last years in order to impose 
such bounds in hyperbolic models. See [9], [1], [6] for the first results of this topic 
and [5] for a numerical version of this kind of problem. 

That is why recently, Berthelin, Degond, Delitala and Rascle ^ proposed a new 
second-order model, which aim is to allow to preserve the density constraint n < n* 
at any time. The main ideas are: 

• modifying the Aw-Rascle model, changing the velocity offset into 

p(n) = (- - , n < n*, 

\n n J 

thus p{n) is increasing and tends to infinity when n ^ n*; 

• rescaling this modified Aw-Rascle model (changing p{n) into ep(ne)) and tak- 
ing the formal limit when e — )• 0"*" . 

This process leads to a limit system on (n, u) which corresponds to the Pressureless 
Gas Dynamics system: 

( dtn + dx{nu) = 0, 
\ dt{nu) + dxinu^) = 0, 

in areas where n < n* . But a new quantity appears, due to the singularity of the 
velocity offset in n = n*. In fact, denoting by p{t, x) the formal limit of ep{n^){t, x) 
when e — 7- C^, we may have p non zero and finite at a point (t, x) such that n{t, x) = 
n* . Thus, the function p turns out to be a Lagrangian multiplier of the constraint 
n < n* . Finally, we obtain the Constrained Pressureless Gas Dynamics (designed 
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as CPGD) model: 

dtn + dx{nu) = 0, 
dt{n{u +p)) + dx{nu{u +p)) = 0, 
< n < n* , p >0 , (n* — n)p = 0. 

The term p represents the speed capabiUty which is not used if the road is blocked 
and that the cars in front imposes a speed smaller than that desired. We refer to [7] 
for more details on the derivation of the CPGD system in the case of the maximal 
density n* being constant. The case where n* depends on the velocity (n* = n*{u)) 
is more realistic (taking into account the fact that the maximal number of cars is 
smaller as the velocity is great, for safety reasons) and is treated in |8]. In [13], a 
numerical treatment of traffic jam is done. 

In this paper, we propose another type of improvement based on the following 
idea: the idea is to use the maximal constraint to model the number of lanes. The 
constraint n* will depend on the number of lanes in the portion of the road. Indeed, 
in a two-lane portion, n* can be twice greater than it is in a one-lane portion of the 
road. This idea simplifies the model dramatically and we no longer need to consider 
as many equations of lanes which makes the modeling and use much simpler while 
reporting the same phenomenon. 

The paper is organized as follows: in the next section, we make a modification of 
the CPGD system to model traffic jams in multi-lane. In section 3, we present the 
dynamics of jams. By approximation with such data, it is used in section 4 to prove 
the existence of weak solutions for any initial data. 

2 The ML-CPGD model 

We consider a piecewise constant maximal density of vehicles, given by 

M 
j=0 

where 

n* € {1,2}, {i^j)i<j<M an increasing sequence of real numbers, 

ro = -OO, TAf+l = +00. 

It means that we set on a road with one or two lanes, the road transitions (change 
of number of lanes) being at points (?'j)i<j<A/- On a one-lane section, the maximal 
density is one (in view of simplification), whereas on a two- lane section, the maxi- 
mal allowed density is two. It is the first improvement of our model: the constraint 
density changes with x to model the fact that there is one or two lanes. Evolu- 
tion equations are given by the Multi-lane Constrained Pressureless Gas Dynamics 
system (designed by ML-CPGD), whose conservative form is 

dtn + dx{nu)={), (2.1) 
dt{n{u + p)/a) + dx{nu{u + p)I^) = , (2.2) 
0<n<n*(x), 'u>0, p>0, (n*(x) - n)p = , (2.3) 
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where the function = /q,(x) is defined by 

I (x)-l ^ if n*{x) = 1, 
\ 1/a if n*(x) = 2. 

The number a > 1 stands for the rate between two-lane velocities and one-lane 
velocities. Thus a single car (we mean a car not into a jam) with speed u on a 
one-lane road will pass to the speed au on a two-lane road. This represents the fact 
that on a two-lane section, the average velocity is higher than on a one-lane (on 
a highway, you drive faster than on a road even if you are alone). The preferred 
velocity depends on the road width according to a. It can also be understood as the 
speed limitation on the various kind of roads. This is the second improvement of our 
model. It only act on the second equation since it is the momentum quantity which 
has to be changed and not the conservation of the number of cars (first equation) . 

Of course, this model can be extended to case with three-lane, four-lane por- 
tion... In the case of three lanes, n* G {1,2,3} and la is replaced by Iot,p{x) = 
1 if n*{x) = l, 

1/a if n*{x) = 2, with /3 > a > 1, a being the rapport of speed between one 
1//3 if n*(rc) = 3, 
and two lanes and /3/a the rapport between three and two lanes. 



3 Clusters dynamics 



In this section, we present some particular solutions (n, u^p) of (2.1 )-(2.3 ) which are 
clusters solutions. For these functions, n = n{t,x) take as only values and n*{x). 
In some sense, they are an extension of sticky particles of |10[ playing a crucial 
role in the proof of existence of solutions for constraint models. They have been 
introduced in [9] and used with various dynamics in fH [6l [71 [!]• 

Let us consider the density n(t, x), the flux n(t, x)u{t, x) and the pressure n{t, x)p{t, x) 
given respectively by 

N 



n{t,x) =n*{x)J2Kit)<x<b,(t), (3.1) 

N 

n{t, x)u{t, x) = n*{x) ^ Ui\^(^t)<x<bi{t), (3-2) 

i=l 

TV 

n{t,x)p{t,x) = n*{x)^pila,{t)<x<b,(t), (3-3) 



i=l 



with 



iVeN*, Ui>0, Pi>0, 
as long as there is no collision and no change of n*(x). That is to say 

ai{t) < bi{t) < a2{t) < b2{t) <...< aN{t) < bN{t) 
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and the number of blocks A'^ is constant until there is a shock or a change of width 
(thus we have N = N{t)). 

This type of picccwisc constant solution writes as a superposition of blocks with 
{n,u,p) = {n*,Ui,pi) constant. Each block evolves according to the interactions 
with the other blocks and the changes of width. 
We have to explain three dynamics: 

• What happens when two blocks collide ? (how to describe a shock) 

• What happens when the road narrows (n*(x) was 2 and becomes 1) ? 

• What happens when the road widens {n*{x) was 1 and becomes 2) ? 

First, let us present some technical properties that will be used in the various cases. 

Lemma 3.1 Let be s,a & [0, +oo[, ip G VQO, +oo[xM), and a, b G C^{]inf{s, a), sup{s, a)[). 
We set 

i-a Mt) 

J{s,a,a,b,u) := / / {dtip{t,x) + u{t)dxip{t,x))dxdt. 

Js J a(t) 

Then we get 

Ma) Ms) 

J{s,a,a,b,u) = / ip{a,x) dx — / ip{s,x)dx (3-4) 

J a{a) J a{s) 

+ 



Proof : We have 
di 







/ (p{t,x)dx 


-J 


Ja{t) 


J a 



I (p{t,b{t)){u{t) -b'{t))dt+ / (p{t,a{t)){a'{t) -u{t))dt. 

J S J S 



fh(t) 

dt(f{t, x) dx + ip{t, b{t)) b'{t) - ip{t, a{t)) a{t), 



thus 

l-a Mt) Ma) Ms) 

/ / dt<p{t,x) dx dt = I ip{a,x)dx— / (p{s,x)dx 

Js J a{t) J a{a) J a{s) 

pa pa 

- I (p{t,b{t))b\t)dt+ / ip{t,a{t))a{t)dt. 

J s J s 



pa po\t) pa pa 

Moreover / / dx<f{t,x)dxdt= / ip(t,b{t)) dt— </?(t, a(i)) and the result 

J s J a{t) J s J s 

follows. □ 



Remark 3.2 We notice that 

J{a, s, a, b, u) = —J{s, a, a, b, u), (3-5) 

J(s, a, b, a, u) = —J{s, a, a, b, u). (3-6) 
// we have a' = b' = u, then 

fb{a) rb{s) 

ip{a,x)dx— I 

i{a) Ja{s) 
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Ma) Ms) 

J{s,a,a,b,u) = / ip{a,x)dx— / ip{s,x)dx. (3-7) 

J a(a) Ja(s) 



If we have a' = u and c is constant, then 



J{s,a,a,c,u) = / (p{a,x)dx— / (f{s,x)dx+ / (f{t,c)u{t)dt. (3.8) 

J a((T) J a(s) J s 

Lemma 3.3 We have the following formulas: 
If a' = b' = c' = u, then 

Ms) 

J{s,a,a,b,u) + J{a,T,a,c,u) = — / ip{s,x)dx 

t\ , , (3-9) 

b{a) pc(T) 



+ / (p{a,x)dx+ / (f{T,x)dx. 

J c((t) J a{T) 

If a' = b' = u and c = b{a) = a{T), then 

J{s,a,a,b,u) + J{a,T,a,c,u) = — / ip{s,x)dx+ / u{t)(p{t, c)dt. (3.10) 

J a(s) J a 

Proof : We have 

J{s, a, a, 6, u) + J(cr, r, a, c, u) 
rbia) Ms) 

ip{a,x)dx— / (p{s,x)dx 

I a{a) J a(s) 



+ 



/a pa 
ip{t,b{t)){u{t) -b'{t))dt+ / ip{t,a{t)) {a{t) - u{t))dt 
J s 

rcir) rcia) 

+ / (p{t,x) dx — / ip{a,x)dx 

J ah-) J a{a) 

+ / Lp{t,c{t)) {u{t) - c'{t))dt+ ip{t,a{t)){a'{t) - u{t))dt 



f6(s) rr 

(p{s,x)dx+ / (p{t,c{t)){u{t) - c'{t))dt 

I a{s) J a 

Ma) /.c(t) 

+ / ip{a, x) dx + / ip{T, x) dx 

J c(ct) J a(T) 

+ / ip{t,b{t)) {u{t) - b' {t)) dt + I ip{t,a{t)){d{t) -u{t))dt. 



Since a' = b' = u, the two last terms vanish and we have 

Ms) 

J{s,a,a,b,u) + J{a,T,a,c,u) = — ip{s,x)dx 



+ / ipit,c{t)) {u{t) -c{t))dt 



Ma) Mt) 

+ / ip{a,x)dx+ I ip{T,x)dx. 

J c(a) J a(r) 



The formulas (3.9) and (3.10) follow. □ 
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3.1 About uniqueness of the dynamics 



In order to work with the most reahstic solution, it is necessary to impose a certain 
number of criteria on the dynamics in question. This discussion also improve the 
paper [7]. 

A single block for which u + p stays constant is a solution, for example the 
function corresponding to the following figure: 



t n 




X 



Remark 3.4 To understand the meaning of the dynamics, for every figure, the 
term {n,u,p) on a zone corresponds to the constant values of the functions on a 
block. 



In fact, in an open subset c]0, +oo[txMa;, where n* is constant, it is very easy 



to see that the dynamic displayed on figure satisfies (2.1)-(2.3), for any value of 
< < Uj. 

Now, remember that the term p represents the speed capability which is not used 
if the road is blocked and that the cars in front imposes a speed smaller than that 
desired. The term is if the density is not n*{x) since in this case the car can go 
to its preferred velocity. Thus there is no reason for a single car to have a nonzero 
pressure term if there is no one before him. And the relation (n*(x) — n)p = do 
not impose p = for the first car of the jam. This is why we assume that the blocks 
satisfy the additional constraint: 

(n*(x) -n(x+))p = 0, (3.11) 

in zones where n* is constant. In this property, we denote by n(x^) the limit, if it 
exists, of n{y) when y — )• x with y > x. With this condition, the dynamics of the 
previous figure is a solution only if Uj = Ui. 

The interpretation is the following : if the first car of the jam has the opportunity 
to use its preferred velocity, it uses it and p becomes zero. If not, p is not necessarily 
zero. 



This is why for various blocks sticking one after the other, the constraint (3.11) on 
p gives the two situations of the above figures. 
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(W2,P2) 



(^^3,0) 





Xq 



(U2,P2) 



(W3,P3) 







Xo 



In fact, an other criteria than (3.11 ) related to the minimization of p > can be 
used. It is clear than choosing p = minimize the p term in the previous described 



situations. When the constraint (3.11 ) cannot be imposed, uniqueness criteria which 



is natural is the minimization of p > 0. 

We now detail the various cases that can appear in the dynamic of clusters. 



3.2 Collision between two blocks without change of width 

In a zone where n*{x) = n* is constant, we consider two blocks {n*,ui,0) and 
{n*,Ur,0), with ui > Ur- Thus, at a time t* > 0, the left block reaches the right 
one, and collide with it. The dynamic is displayed in the following figure. 



t n 




X 
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The density n{t,x), the flux n{t,x)u{t,x) and the function p{t,x) are locahy given 
respectively by 

= / ^*K{t)<x<k(t) + n*tarit)<x<br{t) if t < t*, 
\ ^*\it)<x<k{t) + ^*K{t)<x<bAt) if 



n{t, x)u{t, x) 

and 



n-*Uilai{t)<x<bi(t) + n*Urla,.{t)<x<br(t) if * < t* , 
'^*'^r\^t)<x<k{t) + n*UrK{t)<x<br{t) if t > t*, 



n{t, x)p{t, x) 



if t < t* 

-n,)]la^(i)<^<fe^(i) if t>t* 



with the linear functions ai, bi, Ur, br, di, bi are given by 

j^ai{t) = jbi{t) = uu ai{t*) = a\ bi{t*) = x\ 

Ji'^rit) = -^brit) = Ur, ar{t*)=X*, br{t*) = b* , 

^di{t) = -^bi{t) = Ur, di{t*) = a*, bi{t*) = x*, 



and 



Ul > Ur- 



The left block obtains the velocity of the one being immediately on its right when 
they collide. We extend this when more than two blocks collide at a time t* , by 
forming a new block with the velocity of the block on the right of the group. 



Lemma 3.5 The previous dynamic satisfies (2.1)-(2.3). 

Proof : Let be an open neighborhood of the shock zone (displayed in the previous 
figure). Then, we have, for any continuous function 5 and any test function ip S 

v{n), 

{ dt{nS{u,p,Ia)) + dx{nuS{u,p,Ia)),V' ) 

f +00 



n{t, x)S{u{t, x),p{t, x), Ia{x)){dt(p + udxip)dxdt 
-n*S{ui,0, /a) J(0, t*,ai, bum) 

—n*S{Ur, Ul — Ur, Iq) oo, di, bl,Ur) 

—n*S{Ur, 0, /q,) J(0, OO, ar, br,Ur) 

px* 

{—n*S{ui,0,Ia) + 'n*S{Ur,Ui — Ur,Ia)) / (p{t*,x)dx. 

J a* 
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For S{u,p,Ia) = 1, we get 

( dtn + dx{nu),(p ) = 0. 
For S{u,p,Ia) = {u + p)Ia, we get 

{dt{n{u + p)Ia) + dx{nu{u + p)Ia),ip) = 0. □ 



3.3 Narrowing of the road without colHsion 

Let us move to the situation where the road narrows {n*{x) was 2 and becomes 1). 
Here, we describe the evolution of a block which undergoes this narrowing. The 
speed will be divided by a. 

The dynamic of the block is exhibited in the following figure. 




The density n{t, x), the flux n{t, x)u{t, x) and the functional p{t, x) are locally given 
respectively by 



(21 



n{t, x) = < 



ai{t)<x<bi{t) 



if t<t*, 



21, 
1 



aint{t)<X<X* + K*<x<bi„t(t) t* <t < t**, 



< di{t)<x<bi{t) 



n{t, x)u{t, x) = < 



if t>t**, 

if t<t*, 



and 



n{t, x)p{t, x) 



'^Uintlai„t{t)<x<x* +'"'int'^x*<x<bint(t) t* < t < t** , 

^'^i\{t)<x<bi{t) li t>t**, 



'^PiKi{t)<x<bi{t) if t < t*, 

'^Pint\i„^{t)<x<x* +PintK*<x<bir,t{t) t* <t< t**, 

.Pi\i{t)<x<bi{t) if t > t**, 
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with 



^ai(t) = ^bi{t) = Ui, ai{t*) = a*, bi{t*) = x* 



d 
dt 
d_ 
dt 



aintif) = Uint, aint{t*) = a*, aint{t**) = X*, 
bint{t) = Uint, bintit*) = X* , bint{t**) = b* , 



^~a,{t) = ^h{t) = Ui, ai{t**) = x*, bi{t**) = b* 
at at 



Lemma 3.6 The previous dynamic satisfies (2.1)-(2.3) if and only if 

, - s Ui+Pi \ ^ Ui+Pi 

Pint = Ui+Pi- Uint, [Uint, Pint) = ^Uint, 2Uint ) , Pi 



a 



a 



with 



Ui + Pi Ui + Pi 

< Uint < — , <Ui < 



2a a 
Proof : We have, for any continuous function S and any function 99 G T^i^), 

{ dt{nS{u,p,Ia) + dx{nuS{u,p,Ia)),V ) 

f +00 

n{t, x)S{u{t, x),p{t, x), Ia{x)){dtip + udx^)dxdt 

-2S{ui,pi, l/a)J{0,t*,ai,bi,Ui) 
-2S{Uint,Pint, l/a)J{t*,t** 

, Oiint , X , Uint) 

-S{Uint,PintJ l)J{t*,t**,X*,bint,Uint) " S{Ui,Pi, I) J{t** , CO , di,bi, Ui 



-2S{ui,pi,l/a) / if{t*,x)dx 



-2S{Uint,Pint, 1/a) 



(p{t* , x)dx + Uint / (p{t,x*)dt 



b* 



S{uint,Pint,'^) \ / (fit** , x)dx - Uint ip{t,X*)dt 



+S{ui,pi,l) / ip{t**,x)dx, 

J X* 



thus we have 

dt{nS{u,pJa) + dx{nuS{u,p,Ia)) 

= 2 {S{Uint,Pint, '^/Ol) - S{Ui,Pi, 1/q)) _^..] - t*) 

+ {S{Ui,Pi, 1) - S{Uint,Pint, 1)) '^[x*,b*]{x)S{t - t**) 

+ {UintS{uint,Pint, 1) " '^Uint S {Uint , Pint , l/o)) T^lt* ,t"]{t)^{x - X*). 
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For S{u,p,Ia) = 1, we get 

dtn + d^{nu) = {uint - 2uint) l[t*^t"]{t)Six - x*). 
For S{u,p,Ia) = {u + p)Ia, we get 

dt{n{u + p)Ia) + dx{nu{u + p)Ia) 

= f {Uint+Pint - Ui -Pi) t[a*,x*]{x)6{t - t*) 
+ {Ui+Pi - Uint -Pint) '^[x*,b*]{x)S{t - t**) 

+ {uint{Uint+Pint) " '^^int ( "'"^+^'"0 ) ,t"] " X*). 



Therefore, {n,u,p) is a solution of (2.1 )-( [2^ if and only if 



Uint — '^Uint 

Uint + Pint = Ui+Pi 

Ui+Pi = Uint + Pint 

[ Uint+Pint = ^(uint+Pint) 



Uint — '^Uint 

Uint + Pint = Ui+Pi 

Uint+Pint = U'^i+Pi) 

[ Ui+Pi = ^{Ui+Pi) 



Since Uint,Pint,Uint and p^^^ are nonnegative, it concludes the proof of lemma. □ 
Now, we can find the dynamics governing a single block (n* , Ui , 0) which under- 



goes a narrowing of the road: according to subsection 3.1 and additional constraint 

Pi = 0, pi = 0, 



(3.11 ), we have 



Pint 



in the relations of lemma 3.6 which leads to 

Pint — Ui Uint J Uint — '^Uintj Uint 



0, 



Ui_ 

2a 



Ui 



Ui 



a 



Finally, the only dynamics compatible with (3.11) for a narrowing is: 

t ' 
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3.4 Enlargement of the road without coUision 

Now we explain what happens for a block when the road widens {n*{x) was 1 and 
becomes 2). 

In fact, the block (which comes with n = n* = 1) becomes a block with n = n* = 
2, but its speed will be multiplied by the parameter a. The dynamic is exhibited 
hereafter: 




The density n{t, x), the flux n{t, x)u{t, x) and the functional p{t, x) are locally given 
respectively by 



n{t, x) 



' '^ai{t)<x<bi{t) if t <t*, 

X*<X<bint(t) if t* <t< t**, 



\i„t{t)<x<x* + 21^* 



911 

'^^di{t)<X<bi{t) 

Ui\i{t)<x<bi{t) 



n{t, x)u{t, x) = < 



if t>t**, 

if t < t*, 



Uint\i„t{t)<x<x'- + '^''^int'^x*<x<bi„t(t) t* < t < t** , 



ai{t)<x<bi{t) 



and 



n{t, x)p{t, x) = < 



if t>t**, 
if t<t*, 



Pi^a^{t)<x<bi{t) 

Pint\^„t{t)<x<x* + '^PintK*<x<bint{t) t* < t < t** , 

l^P^\{t)<x<k{t) if t>t**, 



with 



^ai(t) = ^bi{t) = Ui, ai{t*) = a*, bi(t*) = x* , 



dt 



a-intit) = Uint, aint{t*) = a*, aint{t**) = X* 



di 



bintit) = Uint, bintit*) = X* , bint{t**) = b* , 
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Lemma 3.7 The previous dynamic satisfies (2.1)-(2.3) if and only if 



Pint = Ui+Pi-Uint, iUint,Pint) = + Pi) " i Pi = 0!{Ui+Pi 

with 

< Uint <Ui+pi, <Ui < a{ui +Pi). 
Proof : We have, for any continuous function S and any function ip £ 'D{Q), 

{ dt{nS{u,p,Ia) + dxinuS{u,p,Ia)),^ ) 

f + OO 

n{t, x)S{u{t, x),p{t, x), Ia{x)){dt(p + udx(p)dxdt 





= -S{ui,Pi,l)J{0,t*,ai,bi,Ui) - S {uint , Pint ,l)J{t*,t**, aint ,X*, Uint ) 
-2S{uint,PintA/a)Jit*,t**,X*,bint, Uint ) 

+2S{ui,pi, l/a)J{t**,oo,di,bi,Ui) 

r-x* 

= -S{ui,pi,l) / ip{t*,x)dx 

/ X* t" \ 

-S{Uint,PintA) i- j (fit* , x)dx + Uint ip{t,X*)dtj 

-2S{Uint,Pint^^/(^) (y^l ^{t** ,x)dx -Uint Lp{t,X*)dt^ 

+2S{ui,pi,l/a) / ip{t**,x)dx, 

J X* 

thus we have, in 2?' ($7), 

dt{nS{u,p,Ia) + dxinuS{u,p,Ia)) 

= {S{Uint,Pint, 1) - S{Ui,Pi, 1)) t[a*,x*]ix)6{t - t*) 

+2 {S{u^,pi, 1/a) - S{uint,Pint, " 

+ (2UintS{Uint,Pint^ 1/") " Uint S (Uint , Pint , 1)) - X*). 

For S{u,p,Ia) = 1, we get 

5tn + (?a.(nii) = {2Uint - Uint) t[t*^t**]{t)Six - X*). 
For S{u,p,Ia) = {u + p)Ia, we get 

dt{n{u + p)Ia) + dx{nu{u + p)/^) 

= (-Umt +Pmi - -pi) l[a*,x*]{x)S{t - t*) 

+ f {Ui+Pi - Uint -Pint) '^[x*,b*]{x)S{t - t**) 

+ (2u^nt ( "-*+^-* ) - Ui„t(ni„t l[t,^t,,]{t)6{x - X*). 
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Therefore, such a function (n^u^p) is a solution of (2.1)-(2.2) if and only if 

Uint + Pint = Ui+Pi 
Ui+Pi = Uint + Pint 

, Uint + Pint = Ol{uint + Pint) 

and we conclude as in lemma |3]6) n 



int 

Uint ~\~ Pint 
Uint + Pint 
L Ui+Pi 



Uint 
Ui +Pi 

a{ui +pi) 
a{ui +pi) 



Now, if a single block {n* , Ui, 0) undergoes a enlargement of the road, we have 

Pi = 0, Pi = 



in the relations of lemma 3.7 which leads to 

Pint = Ui Uinti {uintiPint, 

with 



f Uint 

V 2 ■ 



aui 



Uint \ 



Ui 



aui 



2aui and then 



< Uint < Ui. 

In this case, we can't impose = since it would imply Uint 
Pint < which is impossible. Then, we use the second criteria of section 3.1 which 
is the minimization of pint in this case. Here pint = is possible, that is the choice 
Uint = Ui, and the dynamics for a enlargement is the following: 



(2,o«,,0) 




This situation is the only case where (3.11 ) can't be additionally asked. The physical 



explanation is that the increasing of the speed of the car going from Ui to aui is not 
instantaneous and has to be in two steps. Thus, in the intermediate state, the car 
is not yet at its preferred velocity and there is still a p term. 

Remark 3.8 We notice that the dynamics for enlargement is exactly the reverse 
process of the narrowing. 



3.5 Compatibility of the dynamics 

Since the previous dynamics are not instantaneous, they can interact before they 
are completed. In this subsection, we present the various compatibilities between 
these dynamics. Note that it is not just a superposition of various cases. In order to 
simplify the presentation, we only show figures that describe the various interactions. 



16 



3.5.1 A train of blocks undergoes an narrowing 




3.5.2 Two blocks collide just before the road ncirrows 

Case with Ui-i > uf 



Case with 
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3.5.3 A train of blocks undergoes an enlargement 



tl 


(2.f .(a-5)"< + «P<) 


(2, a{ui + p<);_o)___^--:::r:^ 






( (2,f,(a-i)ui) 













3.5.4 Two blocks collide just after the road widens 




Here, we have ^ > tti+i, thus 

aui - Uj+i > Uj - 2ui+i > 0. 
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3.5.5 The road follows 1 — t- 2 — )• 1 faster than the block 



to 



1 1 






(2,f ,(a-i)«,,) 












3.5.6 The road follows 2 1 2 faster than the block 



to 



1 1 










(2,|j,(l-^)«0 / 
























3.6 Block solutions and bounds 

Using the above sections, we are able to state some results on the dynamics of 
blocks. 

Remark 3.9 The velocity u is assumed to be extended linearly in the vacuum ( areas 
such that n = between two successive blocks. Moreover, we assume that u is 
constant at ±00. But concerning p, the constraint {n* — n)p = implies that p = 
in the vacuum, and at ±00. Thus, the computations of total variation in x of u and 
p are different. 

The previous computations show the following results: 

Theorem 3.10 With the various above dynamics, the quantities n{t,x), u{t,x) and 
p{t,x) defined by (f!Op-p3[) and Remark\3^are solutions to ( [2?T| ), ([23]). 

We can also establish some bounds on these solutions: 
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Proposition 3.11 We still denote by n{t,x), u{t,x) and p{t,x) the functions of 



(3.1}-(3.3) and Remark 3.9. These functions satisfy the maximum principle 

< u{t,x) < 2a [esssupyU^{y) + esssupyp^{y)) , (3-12) 

< p{t,x) < 2a {esssupyU^{y) + esssupyp^{y)) . (3.13) 

If we assume furthermore that the initial data in the blocks and p^ are BV 
functions, then we have, for all t G [0,T], 

TVK{u{t, .)) < AaM (ry^(u°) + Ty^(p°) + |k°||L-) , (3.14) 

TVK{p{t, .)) < 4aM (rF^(nO) + Ty^.(/) + ||n°|Uoo) , (3.15) 
for any compact K = [a, b] and with 

K = [a — t (esssupyU^), b — t (essinfyU^)], 

where TVk (resp. TV^) denotes the total variation on the set K (resp. K), and 
M is the number of road transitions (supposed to be finite). 

Proof : We treat some examples which represent the critical cases. In these cases, 
we compute the total variation on M to simplify the presentation. 

• Case of collisions without change of width: 

We obtain the bounds corresponding to the classical CPGD model (like in [8|). 

We assume the following dynamics: at time t = 0, there are N blocks (denoted 
by -Bi, . . . , -Bat) with velocities u\> u^> ■ ■ ■ > (which is the case with the 
most collisions) and pressures p\, . . . > 0, thus 

N-l N 

TV{u') = Y^\u^-uU, TV{p') = 2Y^pl 

i=l i=l 

Let t > such that in the time interval [0,t], the j first blocks Bi, . . . , Bj 
collide successively at ti < • • • < tj-i < t for instance (i.e. Bi collide with 
Bi^i at the time ti for all 1 < i < j — 1) and then the q — j + 1 following blocks 
Bj, . . . , Bq collide at the same time tj, with tj„i < tj < t. At the time t, the 
last N — q + 1 blocks Bq, . . . , Bj\f have not collided yet. 

We have the relations 

yk€{l,...,j}, ViG{l,...,iV}, pI^ +ul'' =p^ + ul 
Vi G {1, . . . , q}, = u^q, = + - u% 



ViG {<?+!,... ,iV}, u^=v^, p^=pP 



Then we get 
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TV{u{t,.)) = \u\' - u^\ + --- + \u/_i-u 



■3 I 



3 



+\ tj tj I . . \ tj tj\ 

\Uj - Uj_^-^ \ H 1- - Uq I 



+ \Uq - Uq+l\ + ■ ■ ■ + \U%_^ - U%\ 



+ \Ug ■i^g+il + ■ ■ ■ + li^Ar-i ^atI 



Now, for p, we have 



TV{p{t, .)) = + \p\' -P^\ + --- + - p7 I 

+ \P^j - Pj+l I H ^ - P« I 

+ p'i + 2(f^^, + ...+p%) 

< ti5-u|) + |u?-«^| + --- + |tiO_i-M°| 
+ 2{p'i+p^^ + ---+p%) 

< 2TV{u^) + TV{p^). 

• Case of enlargement of the road without collision: 

We assume the following dynamics: at time t = 0, we consider two blocks 
Bi = {ui,Pi) and B2 = (1*2,^2)' ^ section of road where n* = 1. 
We have 

TV{u^) = \4 - u% TV{p') = 2(p0 

At time ti > 0, the block Bi reach the two-lane section, and undergoes the 
change of width during the time interval [ti,t2]- Then, later in the interval 
[^35*4] (with ts > t2) the block B2 enter in the two- lane section. 
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For all t €z]ti, t2[, we have (with the notations of section 3.4) 

TV{u{t,.)) = \u^ - Ui^intl + \ui^int - UlMtl 

= |n^-ui,i„t|H — 

^ I Ol I I I I '^'^,int 

< - nil + liii - ui,int| H 

Since ui^int < we obtain 

TV{u{t,.)) < |nO-n?|+n?-^ < TV{u'') + \\u^l^ . 
Moreover 



TV{p{t, .)) = 2p^ + + \pi^int - Pl,int\ + Pl,int, 



but 



Pl,int -Pl,int = (a - +Pi) + — > 0, 

thus 

TVipit,.)) = 2ipUpi:^t) = 2(pO + aK+p?)-^), 
and we deduce 

TV{p{t,.)) < 2a{p1 + p'^) + 2au1 < aTV{p^) + 2a\\u^\\Loo. 

For all t £]t2, t^l, we have 

TV{u{t,.)) = |u^-au?| < |n^-u?| + |(l-a)n?| 

< TV{u^) + ia-l)\\u°\\Loo, 

and 

TF(p(t, .)) = 2pl + 2ap\ < aTV{p^). 
For all t G]t3, t4[, we have 

TV{u{t,.)) = \u2,int - U2,int\ + \u2,int " au\\ 
U2,int , ,U2,int Qi 

= -^ + \—-»u,\ 
Since U2 int < n^, we obtain 
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rn^rt \\ ^ U2 int U2 int , i Oi 

TV(u{t,.)) < — ^ ^ \- + \au2 — auil 

< ary(u°) + a||u°||Loo. 

Moreover 

TV{p{t, .)) = P2,int + \P2,int " P2Mi\ + P^Mt + ^OPl 

= 2(a(n0+p0)_!^+«p0)^ 

and we deduce 

TV{p{t,.)) < 2a(p0+p0) _^2an^ < aTV{p°) + 2a\\u°\\Loo . 

At least, for t > ^4, we have 

TV{u{t,.)) = \au^-au^\ = aTV{u^), 

and 

TV{p{t,.)) = 2apl + 2apl = aTV{p°). 
Finally, the bound is 

TViu{t, .)) < a{TViu^) + Wu^Wloo), 
TV{p{t, .)) < a{TViu^) + 2||nO||Loo), 

for all t > 0. 

In the general case (if we follow N blocks along the time), we shall obtain 
the same bound because only one block at a time undergoes every change 
n* = 1 ^ 2. 

But it is possible that many blocks undergo this enlargement together at dif- 
ferent places. That is why the general estimate is the following: 

TV{u{t,.)) < a{TV{u°)+M\\u^\\Lo^), 
TV{p{t,.)) < a(ry(pO) + 2M||uO||L.o), 

where M is the number of lane transitions. 
Case of narrowing of the road without collision: 

The computations are similar to the previous case. With the notations of 



section 3.3, we have 
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For all t e]ti,t2[, 

TV{u{t,.)) = - Ui^intl + \ui,int - Ul^\ 

< 1^^° - ^^ll + - Ul,int \ + Ul,int- 

Since ui^int < we obtain 

TV{u{t,.)) < lu^-ull + u'i < TV{u^) + \\u^\\l'>o. 

Moreover, 

TV{p{t, .)) = + Plaint + \pi,int - Pljiul + Pl^t, 

but this time, pY^ < Pi,int, thus 

TV{p{t,.)) = 2(p0+pi,i„i) = 2(p0 + n0+p?-«i,i„t), 

and we deduce 

TVip{t,.)) < 2(p?+p0) + 2u? < ry(/) + 2||nO||Lco. 

For all t E]t2,ts[, we have 

Tv{uit,.)) = < l^^2-^2l + l(i--Kl 

a a 

< TV{u^) + {l--)\\u^L^, 
a 

and 

TV{p{t,.)) = 2p%+'^p\ < TV(jP). 

a 

For all t Gjts, t4[, we have 

TV{u{t, .)) = \u2,int - U2^t\ + - ^u\\ 

= U2,int + \'2u2,int - ^ul\ 
< U2,int + \2u2,int - -ull + I"""" " 

a a a 



Since tt2,mt < U2, we obtain 
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TViuit,.)) < U2,int + 2ul-2u2,int + {2--)4+-\4-4\ 

a a 



a 



Moreover 



2 

TV{p{t, .)) = P2,int + \P2,int " P2,int\ + P2,int + -Pl 
= 2(«0+pO-U2,int + ^P?), 

and we deduce 

TV{p{t,.)) < 2(pO+pO) + 2tiO < TV{p^) + 2\\u^Lo 
At least, for i > ^4, we have 

TV{uit,.)) = = -ry(t.°), 

a a a 



and 



TVipit,.)) = -pU-Pi = -TV{p% 
a a a 



Finally, the bound is 

TV{u{t,.)) < rF(7x°) +4||u°||loc, 
TV{p{t,.))<TV{u'^) + 2\\u'^\\Loo, 

for all t > 0. 

Now I III" 1 1 LOO can appear on every lane transition and the estimate is then 

TV{u{t,.)) < Ty(«0) + 4M||n°||Loc, 
TVip{t,.)) < TV{p^) + 2M\\u^\\lo., 
where M is the number of lane transitions. 

The general situation is a superposition of these cases and it gives the Propo- 
sition. □ 
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4 Existence of weak solutions 



In this section, we prove the existence of weak solutions using previous clusters 
dynamics, an approximation lemma of the initial data by these sticky blocks and a 
compactness result. 

4.1 Approximation of the initial data by sticky blocks 

We first start by proving the following approximation lemma of initial data. 

Lemma 4.1 Let n° G L^iR), u^,p^ G L°°(M) n BV{R) such that < n° < n*{x), 
< u^, < and (n*(x) — 'nP)p^ = 0. Then, there exists a sequence of block initial 
data {n^,UkTPk)k>i such that for all A; G N*, 

[ nl{x)dx< [ n°ix)dx, (4.1) 
essinf < < esssup n°, essinf p^ <p\< esssup p^, (4.2) 

TV{ul) < TV{u^), TV{pI) < ry(/), (4.3) 

and for which the convergences n^^n^, n^u^^n^u^ and n^p^^n^p^ hold in the 
distribution sense. 

Proof : The proof is widely inspired from the ones in [3] and [8], but here, n* is 
piecewise constant, and constant at ±00. 
Up to a negligible set, we can write 

jez 

where Ij =]aj,aj+i[ is a bounded interval, n*{x) = n*^ for x G /j, and G {1,2} 
(the assumption n* constant at ±00 implies that the sequence {n*)ji=z is station- 
nary) . 

For all A; G N*, we can divide (up to a negligible set) each interval Ij like this: 
fc-i 

= U ]4? ' 45+1 [' 4? = + ~ "j') ' « = 0, . . . , /c. 

1=0 

For j £Z, k £ W, and < i < - 1, we set 



J'* n* L,(k) 



Since < n° < n*, we have < rn~j^^ < !!]£!^(M_ thus 
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We set 



k fc-l 



j=—k 1=0 



Obviously satisfies (4.1). 
Moreover, we can notice that 



nl = a.e. on Ja}/,a}^i[ 



= a.e. on Ja}/,a}4i[, 



and 



n * 1 (k) (k) r 

nl = nj a.e. on ]a) J, a) <= 
We also define 

k fc-l 



n * 1 (k) (k) r 

n'' EE n* a.e. on Jaj^j , aj^Z+ii- 



^fcla^j^^fcla;) = V y^^^^yii^c^) „(*) ,™(fe)r(a;), 

j=—k 1=0 



where uj*^^ = essinf u^, which makes sense because S BV{M.). We have 



la*''' a'''' 



a.e.xG]ag\agVi[, nO(x) / 



^fc(^) = 



We extend li^ linearly in the vacuum (areas where = 0) and at infinity, as in 
Remark 13.91 

Thus, areas where = have no influence on the total variation and we have 



TV{ul) 



(fc) _ (fc) I 

^-fe,0 ^-fc,ll 



/fc) _ ^(k) I 



u 



(k) I 
-fc-+l,Ol 



+ l^-fc,fc-l 

I I (fc) _ (fc) I I ,1 (fc) _ (fc)| 

"I" l^-fc+1,0 ^-fc+l,ll "I + l%-l,fc-l %,ol 



I (fc) _ (fc), I (fc) _ (fc) I 

l%,0 ""fc.ll "I ^ l%,fc-2 %,fc-ll 



< 



TV, 



which shows that satisfies (4.3). We also have (4.2). 
For any test function ip G 'D{M), we have 

j nl{x)ip{x)dx = ^^n* / ip{x)dx 

^ |j|<fc i = '^3,i 



fe-1 



m 



(ky 



|j|<fc i=0 



(fc)^ 
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with a'!'} < f^''} < a~!^] + ni'!'} (if m*;'!-^ ^ 0). Thus, we can rewrite 

|i|<fci=0 

Let Jo G N* such that supp{ip) C | | Ij^ (it is possible because inf {meas{Ij)) > 

\j\<jo 

Then we have, for ah k > jo, 



We also have 



n''(a;)(/?(x)(ix 



E E , 

\j\<jo i=0 •'"■j, 



'nP{x)(f{x)dx. 



Thus 



< 



< 



< 



/ 'nP{x)(p{x)dx — / n^{x)(p{x)d, 
Jr Jul 

k-l .qC^) k-1 

E E /;;"'n°(^)i^(x)-^(a55)idx+ ii^'iu E E 



fc-1 * (k)' 



\3\<30 i=0 
fe-l 



k—\ C^) fc— 1 

ii^'iioo E E-.* [T(--^)dx+ ii^'iu E E-5-f 



\j\<30 i=0 

k-l 



\j\<jo i=0 



n^'Woc E E( 

\j\<jo i=0 ^ 



f meas{Ij) 



< C{ip,jo) X -. 
Moreover, we have similarly 

/n2(.)„«(x).(x).x 

Jr 



k-l / „a^k) 



«°(x)«SV(«S')*'+^''£V'kS' 



and 



fc-l /.a'.';'^^ 

nP{x)u^{x)ip{x)dx = / nP{x)vP{x)(f{x)dx . 



\j\<joi=0'"'j.i 
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Thus 



/ n^{x)u^{x)(p{x)dx — / •n!l{x)u^{x)(p{x)da 
Jr Jilt 



fc-1 

£ EE, 

\j\<jo i=0 "-i 



n^{x)uf)\ip{x) - ^{afj)\dx 



k-l 



{ky 



+ ii^'iu E E«;'.'^ 

\j\<30 i=0 



k-l „„(*=) 



+ E E 

\j\<jo i=0-^°-j 



n'{x)\u\x)-u^^^Mx)\dx 



< C(^,io) X ll^^olloo X ^ + 2||(^||oo E 

\j\<jo i=0 



\u\x)-uf}\dx. 



Therefore we just need to show that the last term vanishes when A; ^ oo. This is 
raised because 



fe-i 

EE 

lil<io i=o 



(fc) N 

;;;-i„»(x)-„widx 



\j\<jo i=0-^S- 



sup u 







inf -0 



(fe) (k) , i„W r 



s E 



\j\<jo 



meas{Ij) 
k 



/k-l 



1=0 



< 



E 

\j\<jo 



< TV{v^)xC{jo)x-. 

We estabhshed that < n°,(^ > — >■ < n°,(^ > and < n^u^^ip > ^ < nPu^,ip >. 
Finally, we define p° the same way as u^^: 

k k-l 

nl{x)pl{x)= E"i^S^^]aW,aW+mW[^^)' 

j=—k t=0 
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where p^^^ = essinf p^ . But in the vacuum (areas where = 0) we set = 0. 

Thus, we have p\ = p^^^ on each interval ]af^i , af^i^ii- In fact, there are two cases: 

• If n° = n* a.e. on ]a^j'^i , a^j^ij^ii, then = n* and pi = 

• Else, it exists a non negligible subset io cja^'^^^ aj^.^^-^[ where < n*-, and 
p° = a.e. on w, which implies p^^^j = 0, and = = p^^^J a.e. on Joj'i^ Qj'i+il- 



We easily deduce that pi satisfies properties (4.2) and (4.3). 



For the convergence < nlpl ,ip > ^ < rrp" ,(p>, the proof is exactly the same 
as nlul □ 

Remark 4.2 The sequence {nl,ul,pl) satisfies the constraint: 

(n*(x)-nO)pO = o, VA: > 1, 
because nl{x) G {0, n*} for all a; £ M. 

4.2 Existence result 

Let us recall the ML-CPGD system: 

dtu + d^{nu) = Q , (4.4) 
dt{n{u + p)Ia) + dx{nu{u + = , (4.5) 
0<n<n*(x), ii>0, p>0, (n*(x) - n)p = 0. (4.6) 

We prove now the existence of weak solutions. The idea is first to approximate 
the initial data in the distributional sense by sticky blocks. These special initial 
data give a sequence of solutions. Then we perform a compactness argument on 
this sequence of solutions. Finally, we prove that the obtained limit is a solution 
for the wanted initial data. The regularity of the solutions are 

n G L°°(]0, +oo[t, L°°(M.) n L^(M.)), (4.7) 

Ti,pGL~(]0,+oo[t,L°°(M,.)). (4.8) 

Theorem 4.3 Let {n^,u^,p^) be some initial data such that 

n° G L\R), u°,p^ G L°°(M) n BV{R), 

with < u^, < p^ , < < n*{x) and (n*(x) — n^)p^ = 0. Then there exists 



{n,u,p) with regularities (4.7), (4.8), solution to the system ( |4.4[ ) — (4.6), with initial 
data {n^,u^,p^). The obtained solution also satisfies 

< u{t,x) < 2a [esssupyU^{y) + esssupyp^{y)) , (4-9) 

< p{t,x) < 2a {esssupyU^{y) + esssupyp^{y)^ . (4.10) 
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Proof : Let n^, u^, {k £ N*) he the block initial data associated respectively to 
nP, vP, provided by Lemma 4A_ For all k, the results of section [s] allow us to get 
{nk,Uk,Pk) solutions of (4.4) — (4.6) with initial data {n^,u^,Pk)^ with regularities 
(4.7), (4.8), and which satisfy the bounds 

< Uk{t,x) < 2a {esssupyu\{y) + esssupyp1{y)) , (4-11) 

< pk{t,x) < 2a {esssupyU^iy) + esssupyp^{y)) , (4-12) 

TVKiukit, .)) < AaM {TV^iuD + TV^{pl) + \\uI\\l^) , (4.13) 

TVKiPkit, .)) < 4aM {TVj,{ul) + TV^{pl) + \\uI\\l^) . (4.14) 
Since (n^) is bounded in L°°, then there exists a subsequence such that 

nfe-n in L-(]0,+oo[xM). (4.15) 

the 



4.1 



Thanks to ( |4.11 ), ( |4.12 ) and the bounds on u^,p^ provided by Lemma 
sequence (ufc) and (pk) are bounded in L°°(]0, +oo[xM), then, up to subsequences, 
we have 

Uk^umLZ{]0,+oo[xR), (4.16) 

Pfe^pinL~(]0,+oo[xM). (4.17) 

Next step is now to prove the passage to the limit in the equation. 
First, for the sequence {nk)k>i, we can obtain more compactness using the following 
lemma and the estimate: 

VT>o, y^eV{R^), vt,se[o,r], v^gN*, 



/ 



{nk{t, x) - rikis, x))ip{x)dx 



< n* sup\\u^\\lo 
k>l 



/ \dx^\dx \ \t — s\ 
Jr J 



(4.18) 



which can be obtained by integrating (4.4). 

Lemma 4.4 Let {nk)k(^N* be a bounded sequence in L°°(]0,T[x]R) which satisfies: 
for all ip G D(Ma;), the sequence [J^nk{t, x)ip{x)dx) is uniformly Lipschitz contin- 
uous on [0, T], i.e. 

3C^>0, VA;eN*, Vs,tE[0,r], 



(nfc(t, x) - rikis, x))Lp{x)dx 



< C^\t - s\. 



Then, up to a subsequence, it exists n G L°°(]0,T[x]R) such that rik ^ n in 
C([0,r],L-(M,)), I.e. 



Vr G L^(]I 



sup 

te[o,T] 



{nk{t, x) — n{t, x))T{x)dx 



0. 
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Proof : Let (</?m)m>i be a countable set dense in 'D{M.x) for the L^-norm, which 
exists because of the separabihty of L^{M.x)- We denote 

gk,m{t) ■■= / nk{t,x)(pm{x)dx. 

The sequence {gk,i)k>i is bounded and equicontinuous in C([0,T],M), thus, the 
Ascoh Theorem entails that it exists an extraction (Ti{k) such that 

9<xi(fc),i^^i in Ci[0,T],R). 

The same applies to {gai(k),2)k>i, thus it exists an extraction a2 such that 

g.,{a,{k)),2 ^ h in C([0,T],M). 

A simple recursion shows that we can build a sequence of extractions am such that 

fl'fTi((72{...<T,„(fc))...)),m ^™ in Ci[^,T],M). 

Therefore, setting cr{k) := cri o • • • o a^^k), we have (by diagonal extraction) 

Vm>l, ga(k),m-^lm in C([0,r],M). (4.19) 

A: 

Now, we can identify the limit Im because since {n^fj^^)}^ is bounded in L°°(]0, T[x]R), 
there exists a subsequence (still denoted by the same way) such that ^^(fc) ^ n in 
L2'*(]0,r[xM). Thus, we have, for all m > 1, and for all G P(]0,r[i), 

/ / n„(k){t,x)'il^{t)y^m{x)dxdt ^ / / n{t,x)il}{t)ipm{x)dxdt, 
Jo Jr ^ Jo Jr 



which rewrites 







ga{k),m{t)'4'{t)dt 



n{t, x)(pm{x)dx I il){t)dt. 



Moreover, (4.19) easily implies that 







ga(k),mit)'4'it)dt -> / lm{t)lp{t)dt 



thus lm.{t) = J^'n{t, x)ipm{x)dx, a.e. t £ [0,T], from which we can deduce 



Vm > 1, 



sup 

ie[o,T] 



("<7(fc)(*, 3;) - nit, x))ipm{x)dx 



0. 



Finally, this convergence stays available for all 99 G 2?(M^), because of the inequality 

sup I {n„(^k)-n){t,x)ip{x)dx 
ie[o,T] Jr 



< sup 

te[o,T] 



(?^<7(fc) - n){t,x)iprn{x)d2 



+ C||(/9 - ^mWh^i 
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where 



C := sup(||nA;||2,oo(]o,T[xR)) + II"-IIl°°(]o,t[xR) < +00. 



k>l 



We conclude that it is also true for T £ L^{M.x) by density, using the same inequal- 
ity- □ 

Following of the proof of Theorem 4.3 According to (4.6) and (4.18), the lemma 



4.4 



applies to the sequence (nfc)fc>i, and thus 

nfc ^ n in C([0, T], L^,(M^)), for all T > 0. 



(4.20) 



As the same, we obtain (integrating ( |4.5| )) an estimate similar to (4.18) for the 
sequence {nk{uk + Pk)Ia)k>i-, thus it exists q G L°°(]0, +oo[xM) such that 



nk{uk+Pk)Ic. ^ q in C([0, T], L^,(M,)), for aU T > 0. 



(4.21) 



Now, the key point of the proof is passing to the limit in the products and is treated 
by the following technical lemma: 

Lemma 4.5 Let us assume that {'^k)k&i is a bounded sequence in L°°(]0,T[xM) 
that tends to 7 in L^^(]0, T[xM), and satisfies for any T G T>(Rx), 



ilk - 7) , x)T{x)dx — > 0, 

k 



(4.22) 



either i) a.e. t G]0,r[ or ii) in L^{]0,T[t). 

Let us also assume that (wfc)fcgp} is a bounded sequence in L°°QO,T[xM) that tends 
to bj in L^^(]0, r[xM), and such that for all compact interval K = [a, h], there exists 
C > such that the total variation (in x) of and u over K satisfies 

VA:GN, TVK{uJk{t,.))<C, TVK{uj{t,.))<C. (4.23) 

Then, "fkOJk^loJ in L^^(]0, T[x]R), as k ^ +00. 

Remark 4.6 This is a result of compensated compactness, which uses the compact- 
ness in X for {oJk)k given by (4.23) and the weak compactness in t for {■~fk)k given 
by (4.22) to pass to the weak limit in the product jk^k- 

Proof : We can refer to |4] for a complete proof, even in the case where 



WkeN, rvK(wfc(t,.)) < C(i + 

which is more general. □ 



1. 



TVK{oo{t,.))<C{l + ^), 



End of the proof of Theorem 4.3 The convergence (4.20) allows to apply Lemma 
with jk = iT'k- Moreover, thanks to (4.13) and the BV bounds on provided 



4.5 



by Lemma 4.1, we can set Uk = itfc in Lemma 4.5 (in fact, the sequence Uk{t, .) is 
uniformly bounded in BV with respect to t, and also u{t, .) thanks to the lower 
semi-continuity to the BV norm). Thus, we have 



UkUk^nu in L2'=^(]0, +oo[xM). 



(4.24) 
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The same applies to the sequences {'Jkj'^k) = {^k^Pk) and {'^k-.^k) = i^kiuk + 
Pk)Ia,Uk)- we have 

nfcPfc-npinL~(]0,+oo[xM), (4.25) 

nkiuk+Pk)IaUk^qu in L^^(]0, +oo[xM). (4.26) 
Furthermore, we easily have 

rikiuk + Pk)Ia^n{u + p)Ia in L2'^(]0, +oo[xM), 

thus q = n{u + p)Ia, and 

nkUkiuk + Pk)Ia^nu{u + p)Ia in L^,(]0,+oo[xM). (4.27) 

We deduce that (n, u,p) satisfies (|4.4[), (|4.5[) in V'QO, +oo[xR), and the constraints 



(li6|). 

The last step is to show that {nP ,p^ ,u^) is really the initial data of the problem, 
according to the weak formulation: 

Vv3G C,°°([0,+oo[txM,), 

{ndtf + nudxf) {t,x)dxdt + / n^{x)(p{0,x)dx = 0, 

Jr 

{n{u + p)Iadt(p + nu{u + p)Iadx^p) {t, x)dxdt 
+ / n°(x)(n°(a;) +p°(x))/„(x)v9(0,x)dx = 0. 







It comes easily, because we have, for all > 1: 

V(^E(:7~([0,+oo[txM^), 



oo 







oo 



(ukdiip + nkUkdxf) {t,x)dxdt + / n^(x)(^(0, x)dx = 0, 

Jr 

{nk{uk + pk)Iadt^p + UkUkiuk +pk)Iadx(p) {t,x)dxdt 
+ [ nl{x){ul{x)+pl{x))Iaix)ipiO,x)dx = 0, 



and we can pass to the limit when k — t- +oo because of the convergences -^"nP, 



'k 

n^u^^ vPv? and n^p^^ vPp^ in I?'(M), and the convergences (4.24), (4.26) and 



(4.27) inL-(]0,+oo[xM). □ 



4.3 Compactness result 

To finalize the paper, we set a compactness result which is contained into the proof 
of the previous existence Theorem. 
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Theorem 4.7 Let us consider a sequence of solutions {nk,Uk^pk) with regularity 
(4.7), (4.8), satisfying (4.4) — (4.6), and the following hounds: 



VA; G N, a.e. g]0,+oo[xM, Q < Uk{t,x) < Ca, 

V/c G N, a.e. g]0,+oo[xM, ^ < Pk{t,x) < C^, 

Vi^=[a,&]cM, VfcGN, a.e. t g]0, +oo[, TVK{uk{t, .)) < C^^m,k, 

\/K = [a, h\ CM, VA; G N, a.e. t g]0, +oo[, TVK{pk{t, ■)) < Ca,M,K, 

with Ca (resp. Ca,M,K) some positive constant depending only on a (resp. a, M 
and K). 

Then, up to a subsequence, {nk,Uk,Pk) ^in,u,p) in L'^^QO, +oo[xU) , where {n,u,p) 



is a solution to the system (4.4) — (4.6). This solution {n,u,p) also satisfies 
a.e. {t, x) g]0, +oo[xM, < u{t, x) < Ca, 
a.e. g]0,+oo[xM, 0<p{t,x)<Ca. 
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